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Abstract 



Exact and general results on the electronic states in one dimensional crystals bounded 
at T and r + L - where L — Na, N is a positive integer and a is the potential period - 
are presented. Corresponding to each energy band of the Bloch wave, there are N — 1 
states in the finite crystal and their energies are dependent on the crystal length L but 
not on the crystal boundary r and map the energy band exactly; There is always one 
and only one electronic state corresponding to each band gap of the Bloch wave, whose 
energy is dependent on the crystal boundary location r but not on the crystal length L. 
This state is either a constant energy confined band edge state or a surface state in the 
band gap. A slight change of the boundary location r could change the property and the 
energy of this state dramatically. 



PACS numbers: 73.21.-b, 71.15.-m, 73.61.-r, 71.18.+y 
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Bloch theorem plays a central role in our current understanding on the electronic 
structures in modern solid state physics. However, any real crystal always has a finite 
size and does not have a hypothetical infinite size or the periodic boundaries which 
Bloch theorem is based on|^. The difference between the electronic structure of a real 
crystal of finite size and the electronic structure obtained based on the translational 
invariance becomes more significant as the crystal size decreases. A clear understanding 
of the properties of electronic states in real crystals of finite size has both theoretical and 
practical significant importance. A straightforward way is to obtain exact solutions for 
crystals of finite size. However, to obtain exact solutions of the Schrodinger equation with 
a general periodic potential for a finite crystal with boundaries has been being considered 
as a rather difficult problem: the lack of translational invariance in the crystal of finite 
size is a major obstacle. It is the use of the translation invariance that greatly simplifies 
the mathematics in solving the Schrodinger equation with a periodic potential. Without 
a such simplification the corresponding problem for finite crystals with boundaries could 
become rather difficult mathematically. Thus most previous theoretical investigations on 
this subject were based on approximate and/or numerical approaches and were usually 
on a specific material and/or on a specific model 0. 

Historically, many of our current fundamental understandings on the electronic struc- 
tures of crystals were obtained through the analysis of one dimensional crystals 0, |, 
Among the most well known examples are the Kronig-Penney model0, Kramers' anal- 
ysis of on the band structure of one dimensional infinite crystals]^, the Tamm's surface 
state P], etc. Kramers' analysis on the solutions of the one dimensional Schrodinger 
equation with periodic potential - based on a differential equation theory approach - pro- 
vided a thorough and general understanding on the band structure of one dimensional 
crystals!, 1,1. 

In this work we present exact and general results on the electronic states in one 
dimensional crystals of finite length L = Na - where a is the potential period and 

is a positive integer - in comparison with the band structure of one dimensional 
infinite crystal, using a differential equation theory approach [|I^]. It also shows that the 
obstacle or mathematical difficulty due to the lack of translation invariance in fact can be 
circumvented. Different from Kramers' work, the understanding on the zeros of solutions 
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of one dimensional Schrodinger differential equation with periodic potential |T0[ plays a 
fundamental role in this work. 

The one dimensional Schrodinger differential equation with a periodic potential can 
be written as 

— y"{x) + [v{x) — \]y{x) = 0, — oo < x < + oo (1) 

where 

v{x + a) = v{x). 

We assume (1) is solved, all solutions are known. The eigenvalues are energy bands en{k) 
and the corresponding eigenfunctions are Bloch functions 0„(fc,a;), where n = 0, 1, 
2, 

For a one dimensional crystal with two boundaries at r and r + L, we assume the 
potential inside the crystal is still v{x) as in (1) and all electronic states in the crystal 
are confined in the finite size of the crystal. We are looking for the eigenvalues A and 
eigenfunctions ip{x,A), which are solutions of 

- y"{x) + [v{x) - A]y{x) = 0, t < x < t + L (2) 

and 

y{x) =0, x<Torx>T + L (3) 

where r is a real number, L = Na and is a positive integer. 

By using a differential equation approach we found that all solutions of (2) and (3) 
can be obtained if all solutions of (1) are known (See Appendix A) and the electronic 
states in one-dimensional finite crystals can be classified into two types: 
1. Inside an energy band, there are A^ — 1 electronic states whose energy is dependent 
on the crystal length L and given by: 

fc(A,)=j7r/L, J = 1,2, ....AT -1. (4) 

Correspondingly, there are A^ — 1 eigenfunctions Aj). Each eigenvalue for this case 
is a function of L, the crystal length. But they all do not depend on the location of the 
crystal boundary r or r + L. For simplicity, we call these states as L-dependent states. 
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although only the eigenvalue of a such state is dependent on only L, the wavefunction of 
a such state is dependent on both r and L. 

The eigenvalues k{Kj) map the dispersion relation k{en) of the original equation (1) 
exactly. Many authors noted that there is an approximate correspondence between the 
bulk energy dispersion and the energy levels in crystals of finite size|TT], |T2|, |T3|, from 



Eq. (4) we see that this is in fact an exact correspondence for the electronic states in 
general one dimensional crystals of finite length. Furthermore this exact correspondence 
does not depend on the crystal boundary r. Pedersen and Hemmer [|T^] investigated the 



non-bandedge electronic states in one dimensional finite crystals with a Kronig-Penney 
model and also found that the energy spectrum is independent on the crystal boundary. 
2. For each band gap, there is always one and only one electronic states, whose energy 
A is given by a necessary and sufficient condition: 

V^(r + a,A) = V^(r,A) = 0. (5) 

(5) does not contain the crystal length L, thus the energy A corresponding to a band 
gap is dependent on the boundary of the crystal r, but not on the crystal length. For 
simplicity, we call these solutions as r-dependent states, although only the eigenvalue of a 
such state is dependent on only r, the wavefunction of a such state is dependent on both 
r and L. The energy A of these r-dependent states can be labeled as A^^2m or Kr^2m+i, 
corresponding to the band gap at = vr/a or /c = 0. Kr^2m is in [e2m(7i"/a), £^2m+i(7r/a)] 
and K,2m+i is in [£:2m+i(0), e2m+2(0)]. 

How these r dependent eigenvalues change as the crystal boundary r changes can be 
seen in Fig. 1, using A,- ! as an example. The details can be found in Appendix B. In 
Fig. 1 the zeros of band edge wavefunctions 0i(O, x) and 02(0, x) are shown as circles. 

If r is equal to a zero of a band edge wavefunction 0i(O,x) or 02(0, x), At-_i is equal 
to £i(0) or £2(0). In any one of these cases, we have an electronic state whose energy is 
the corresponding band edge energy of (1), independent on the crystal length L. Zhang 
and Zunger first observed a state with such behavior in their numerical investigations 



on Si quantum films||TT|. Franceschetti and Zunger also observed a such state in their 
calculations on the free standing GaAs quantum film[|l^. The author has pointed out 
that the existence of the band edge states with such behavior in fact could be quite 
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general in one dimensional symmetric finite crystals |T5|. From this work we understand 
that for the electronic states in one dimensional finite crystals, without any requirement 
on the symmetry of the periodic potential or of the system, as long as a boundary is a zero 
of a hand edge wavefunction of (1), there always is a confined state whose energy is the 
band edge energy, independent on the crystal length L. Due to the inversion symmetry 
of the finite crystals treated in for each band gap there is always one band edge 
wavefunction which is zero at the lattice boundaries. 

If r is not a zero of either one of the two band edge wavef unctions, i.e, r is in the 
dotted line region or in the dashed line region in Fig. 1, then the electronic state ip{x, A) 
will have the form of either e^^^^^pi{x, A) or e~^^^^^p2{x, A) with an energy A inside the 
band gap. This is a surface state. 

Many years ago, by using a Kronig- Penney model TammQ] showed that the termina- 
tion of the periodic potential at the surface of a crystal can cause a surface state existing 
in the band gap for each band gap of the Bloch wave. More than sixty years later, Zhang 



and Zunger [pJ]] and Franceschetti and Zunger fL^ observed the existence of the constant 
energy confined band edge state in their numerical calculations. Now we understand that 
in the one dimensional case, the surface state in the gap and the constant energy confined 
state at a band edge are two different results of the termination of the periodic potential, 
depending on whether the boundary r is a zero of a band edge wavefunction. It may be 
noticed that a Tamm's surface state is an extra "added" state, the r-dependent states in 
this work comes from a unified solving of the all states. 

An interesting question is: Why a finite one dimensional crystal of two ends can 
only have at most one surface state in each gap? This is due to that the two ends of 
a finite crystal in general are not equivalent, except in a symmetric finite crystal [0]. 
A symmetric one dimensional finite crystal requires: (1) The crystal potential has an 
inversion symmetry center. (2) The two ends of the crystal are symmetric to the same 
symmetry center and thus are equivalent. If any one of these two is not true, then the 
two ends of the finite crystal are in fact not equivalent. For a symmetric finite one 
dimensional crystal, it has one constant energy confined band edge state, rather than a 
surface state, in each gap. That is the result obtained in |T3 |. 

A little further thinking can also help to understand this unexpected result. For 
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simplicity we assume a monoatomic crystal of L = Na here. If there are — 1 bulk 
states in each band, and there are altogether atoms, there is only one, not two, more 
state which can be the either the surface state or at the band edge. In the case of = 1, 
there is one state for each energy band. Similarly, one monolayer of atoms also can have 
only one two-dimensional energy band for each atomic state, even it has two surfaces. 
The author has also worked out a similar problem on three dimensional finite crystals [|T7 



For a finite crystal of parallelogram shape which has A'^iai, N2aL2 and N^a^, forming the 
sides which meet at a corner, for each energy band there are (A^^i — 1)(A^2 ^ 1)(^3 — 1) 
bulk states, (A^i - l)(A/'2 - 1) + (A/'a - l)(A/'3 - 1) + (A^'s - l)(A^i - 1) surface states, 
{Ni — 1) + (A''2 — 1) + (A's — 1) "side states" and one "conner state" in spite of that the 
crystal has six faces, twelve sides and eight corners: The sum of these numbers is exact 
N1N2N3. Any change of these numbers will not give the correct sum. 

A slight change of the boundary location r can change the properties of the r de- 
pendent state dramatically, this can be clearly seen in Fig. 1: If r is in the region 
corresponding to a dotted line, the surface state is near one end of the crystal. If r is in 
the region corresponding to a dashed line, the surface state is near the other end of the 
crystal. If r is a zero of a band edge wavefunction (either a solid circle or an open circle), 
then r + L is also a zero of the same band edge wavefunction and thus the boundary r 
dependent state is a constant energy confined band edge state. 

As an example, in Fig. 2 is shown a comparison between the the energy bands and 
the eigenvalues of the electronic states in an one dimensional crystal of length L = 8a. 
In Fig. 3. are shown the energies of three electronic states in crystals of finite length 
near or in the lowest band gap of the Bloch wave as functions of the crystal length L. 

In summary, by using a differential equation theory approach, we have obtained exact 
and general results on all electronic states in one dimensional crystals of finite length in 
comparison with the electronic states in a periodic potential obtained from (1). For one 
dimensional crystals bounded at r and r + L, there are two different types of electronic 
states: There are A^— 1 states corresponding to each energy band of (1). Their eigenvalues 
A are given by (4), thus are dependent on the crystal length L but not on the crystal 
boundary location r and map the energy band exactly; There is always one and only one 
electronic state corresponding to each band gap of (1), whose eigenvalue A is dependent 
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on the boundary location r but not on the crystal length L. Such a r-dependent state 
can be either a constant energy confined band edge state (if the crystal boundary r is a 
zero of a band edge wavefunction of (1)) or a surface state in the band gap (if the crystal 
boundary r is not a zero of either band edge wavefunction of (1)). A slight change of the 
crystal boundary location could change the properties and the energy of this boundary 
dependent state dramaticallyfTSl] . 

The results obtained in ||15| are a special case of the more general results obtained 
here. 
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Appendix A 



The Schrodinger differential equation (1) can be considered as a special case of the 
more general Hill's equation |1TU[. This work is mainly based on the general properties of 



solutions of Hill's equation. 

We are mainly interested in the cases for which there is always a gap between two 
consecutive energy bands of (1)||19|. For these Bloch waves, the band edges £^(0) and 



Snij/a) occur in the order 

£0(0) < £o(7r/a) < ei(7r/a) < £i(0) < £2(0) < e2{'K/a) < £:3(7r/a) < £3(0) < £4(0) < 

The band gaps are between e2m(7r/a) and £:2m+i(7r/a) or between e2m+i(0) and < e2m+2(0). 

Suppose yi{x, A) and 2/2(3^, A) are two linearly independent solutions of (1), in general, 
a solution of (2) and (3) if it exists, can be expressed as 

?/'(a;. A) = /(x, A), ifT<x<T + L 

= 0, ifx<Torx>T + L 

where 

/(x. A) = ci?/i(x. A) + 02^/2(2;, A) (A.l) 
is a non-trivial solution of (1) and satisfies 

/(r,A) = /(r + L,A) = 0. (A2) 

The non-trivial solutions of (2) and (3) can be found through (A. 2) based on the general 
properties of solutions of (1). 

The properties of linearly independent solutions of (1) are determined by a real num- 
ber -D(A) called its discriminant : Suppose ?7i(x. A) and 172(2;, A) are two linearly in- 
dependent solutions of (1), which satisfy the initial conditions 771 (0, A) = 1, ?7^(0, A) = 
0; 772(0, A) = 0, 7/2(0, A) = 1. The real number -D(A) defined by -D(A) = rii{a, X)+ri2{a, A) 
is called the discriminant of (1). 

As A increases from —00 to +00, the value of -D(A) changes as the following (m = 0, 
1, 2, ....)[|IU],0,i: 
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(i) In the intervals [£2m(0), £:2m(7r/a)], D{\) decreases from 2 to -2. 

(ii) In the intervals [e2m+i(7i"/a.), £^2m+i(0)], -D(A) increases from -2 to 2. 

(iii) In the intervals (-oo,£:o(0)) and (£:2m+i(0), £2m+2(0)), D{X) > 2. 

(iv) In the intervals (e2m(7r/a), £:2m+i(7r/a.)), -D(A) < —2. 

The existence of and the properties of non-trivial solutions A and f{x,A) in (A. 2) 
can be determined on this basis. 

For solutions of (2) and (3), both the permitted and the forbidden eigenvalue ranges 
of (1) cannot be excluded. We should consider solutions of (A. 2) for A in (— cxo, +oo). 
However, according to the Theorem 3.2.2 of |jlO[, there is not a nontrivial solution of (A. 2) 



for which A in (— oo, £:o(0)]. Thus we need only to consider A in {eo{0), +oo). Depending 



on A, the value of -D(A) may have five different cases flO 
Case A. \D{X)\ < 2. 

In this case A is inside an energy band of (1). Two linearly independent solutions of 



(1) can be expressed as [TO 



y,ix, A) = e^^W^pi(x, A), y^ix, A) = e-^'=W^p2(x, A), 
where k{X) is a real number depending on A and 

< k{X)a < TT 

and pi{x, A) and P2{x, A) have period a: pi{x + a, A) = Pi{x, A). All A;(A) and Pi{x, A) are 
functions of A. 

If there is a non-trivial solution f{x,A) of (A. 2), simple mathematics gives that we 
must have either 

gifc(A)L _ g-ifc{A)L _ (AAl) 

or 

cipi(r,A) = 0, and 02^2(1", A) = 0. {A.A.2) 

We can easily prove that neither pi(r, A) nor P2{t, A) can be zero. Suppose pi(r, A) = 0, 
we must have pi{t + a,A) = and thus ?/i(r, A) = yi(r + a, A) = 0. Then according 
to theorem 3.1.3 of [|lOl, A must be in [£:2m(vr/a), £2m+i(vr/a)] or in [£:2m+i(0), £:2m+2(0)], 
in which |-D(A)| > 2. This is in contradictory to |-D(A)| < 2 for the Case A. Similarly 
P2(t, A) can not be zero. Thus no non-trivial solution obtained from {A.A.2) exists. 
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Therefore non-trivial solutions in this case can only be obtained by (A.A.l). Note 
[A.A.I) does not contain r. The non-trivial solutions can be obtained if 

fc(A)L = JTT, 

Thus in each energy band k{en), there are — 1 values of A^, where j = l,2,....A^ — 1, 
for which 

k{A,) = J 7t/L. 

This is Eq. (4) in the text. 
Case B. D{X) = 2. 

In this case A is at a band edge at A; = 0: A = e2m+i(0) or A = e2m+2(0). 



Two linearly independent solutions of (1) can be expressed as[|TO 



yi{x, A) = pi{x, A), y2{x, A) = xpi{x, A) +P2{x, A), 



and pi{x, A) and P2{x, A) have period a. 



Due to the Sturm Separation Theorem [|20|, the zeros of A) are separated from the 
zeros of P2{x, A). Simple mathematics leads to that the existence of non-trivial solution 
(A. 2) in this case requires 

Pi(r,A) = and cs = 0. (A.B.l) 



Case C. D{X) > 2. 

In this case A is inside a band gap at k = 0: e2m+i(0) < A < £:2m+2(0). Two 
linearly independent solutions of (1) can be expressed as[p!0| 



1/1 (x, A) = e^^Xa:, A), ^^(x, A) = e-^W>2(x, A), 

where /3(A) is a positive real number depending on A and pi{x, A) and P2{x, A) have period 
a. 

Again due to the Sturm Separation Theorem [pO|] , the zeros of pi{x, A) are separated 



from the zeros of p2{x,X). If there is a non-trivial solution f{x,A) in this case, simple 
mathematics from (A.l) and (A. 2) gives that we must have either 

Pi{t,A) = 0, and C2 = (A.C.l) 
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or 

P2(r,A) = 0, and Ci = 0. {A.C.2) 

Note that {A.C.I) and {A.C.2) cannot be true simultaneously. 

From the discussion on Case B and Case C, we can see that if we have a non-trivial 
solution of (A. 2) for a band gap at A; = (e2m+i(0) < A < £:2m+2(0)), either one of 
(A.B.I), (A.C.I) or (A.C.2) must be true. Since all functions Pi{x, X) in (A.B.I), (A.C.I) 
and (A.C.2) are periodic functions, we always have /(r + a. A) = if we have /(r. A) = 0. 
Therefore the following equation is a necessary condition for having a solution A in (A. 2) 
for a band gap at k = 0: 

f{r + a,A) = /(r,A) = 0. (A3) 

It is easy to see that (A. 3) is also a sufficient condition for having a solution (A. 2): From 
(A. 3) one can obtain /(r + n a. A) = 0, where n = 0, ...N. 

Corresponding to a band gap at = vr/a, the Case D (-D(A) = —2) and the Case E 
{D{\) < —2) can be very similarly discussed as the Case B and the Case C. The major 
difference is there we have semi-periodic functions {si{x + a. A) = —Si{x,A.)) instead of 
periodic functions. We are led to the same equation (A. 3) as a necessary and sufficient 
condition for having a solution (A. 2) for a band gap aX k = ti / a. Thus (A. 3) is a necessary 
and sufficient condition for having a solution (A. 2) corresponding to a band gap. (A. 3) 
is essentially the same as Eq. (5) in the text. 

The Theorem 3.1.3 of |]10[ indicates that for an arbitrary real number r, there is always 
one and only one A for which (A. 3) is true for each band gap [e2m{'^ / o) , e2m+i{T^ / o)] or 
[e2m+i(0) , ^2m+2 

Therefore for any real number r there is always one and only one A, which is a 
solution of (A. 2) and dependent on r hut not L, for each hand gap [e2m{T^ / o) , e2m+i{'^ / o)] 

or [£:2m+l(0),£:2m+2(0)]. 
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Appendix B 



How these r-dependent eigenvalues A change as r changes can be obtained from 
differential equation theory. Here we discuss A,- 2m+i as an example. AT-_2m can be very 
similarly discussed. 

According to the theorem 3.1.2 of [0, 02m+i(O,a;) and 02m+2(O,a;) have exact 2m + 2 
zeros in [0,a). Then according to the Sturm Comparison Theorem |pO|, the zeros of 
02m+i(O,x) and 02m+2(O,a;) must be distributed alternatively: There is always one and 
only one zero of 02m+2(O,x) between two consecutive zeros of 02m+i(O,x), and there is 
always one and only one zero of 02m+i(O, x) between two consecutive zeros of 02m+2(O, x). 

We treat r as a variable and let r go continuously from a (any) zero Xi^2m+i oi 
02m+i(O,x) to Xi,2m+2, the zcro of 02m+2(O,a;) that is next to Xi,2m+i- Since A^_2m+i 
is a continuous function of t[|1^], the corresponding AT-^2m+i will also go continuously 
from £:2m+i(0) to £2m+2(0). Similarly if r goes continuously from Xi,2m+2 to X2,2m+i, 
the next zero of (p2m+i{0,x), the corresponding AT-2m+i will also go back continuously 
from £:2m+2(0) to e2m+i(0). In Figure 1 is shown A^-^i as function of r in the interval 
[xi^i, xi^i + a], where is a zero of 0i(O, x). 

If we consider as r goes from Xi,2m+i to Xi,2m+2 and then to X2,2m+i, the corresponding 
AT,2m+i as function of r goes up from £:2m+i(0) to £:2m+2(0) and then back to £:2m+i(0) as 
a basic undulation, in this basic undulation the function f{x,A) has different forms. We 
know for any solution (A. 2) in each band gap [e2m+i(0), £:2m+2(0)], either one of (A.B.I), 
(A.C.I) or (A.C.2) must be true. When A = e2m+i(0) or A = e2m+2(0), (A.B.I) is true. 
For the two sections (end points excluded) of this basic undulation either (A.C.I) or 
(A.C.2) is true. In one section f(x, A) has the form cie^*-"^-* ^ Pi{x, A) (A.C.I is true) and 
in the other section f{x,A) has the form 026"^'-^^ ^ p2{x,A) (A.C.2 is true). In which 
section it has which form is dependent onv{x) in (1). In Fig. 1 the two sections of a basic 
undulation are shown as a dashed line and a dotted line, indicating two different forms 
of /(x. A). Since in the interval [0, a), both 02m+i(O,a;) and 02m+2(O,a;) have exactly 
2m + 2 zeros, then in general AT-2m+i as function of r will always complete 2m + 2 
basic undulations in an interval of length a. Similarly A^- 2m as function of r will always 
complete 2m + 1 basic undulations in an interval of length a. 
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A function with the form of 
is mainly distributed near either one of the two boundaries of the finite crystal, due to the 
exponential factor. Thus these r-dependent states are in fact the surface states introduced 

by the termination of the periodic potential. These surface states are introduced into the 
band gaps when the boundary r is not a zero of either band edge wavefunction of the 
original Bloch waves. Of course, for the band gaps at k = vr/a {e2m{'^/o) < A < 
^2m+i{'^ / o)) 1 a surface state has the form of either cic'^^^^ ^ •Si(a;, A) or C2e~^'^^'^ ^ S2{x, A) 
inside the crystal, where Si{x, A) and 52(2^, A) are semi-period functions. 
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Figure Captions 



Fig. 1. A^ i as function of r in the interval + a]. The zeros of 0i(O, x) are 

shown as sohd circles and the zeros of 02 (0, x) are shown as open circles. Note that A^^i 
completes two basic undulations in [xi^i, Xi^i + a]. The dashed lines and dotted lines can 
be considered as the surface state is located near two different ends of the finite crystal. 

Fig. 2. A comparison between the energy bands en{k) of (1) (solid lines) and the 
eigenvalues A of the electronic states in a crystal of length = 8a (solid circles, L- 
dependent; open circles, r-dependent). Note that: the L-dependent eigenvalues map the 
energy bands exactly and satisfy (4); the r-dependent eigenvalues are in a band gap or 
at a band edge of (1). 

Fig. 3. The eigenvalues of three electronic states in finite crystals near or in the 

lowest band gap of (1) as functions of the crystal length L, while r is fixed. Note the 
eigenvalue of the r-dependent electronic state (open circles) in the gap is independent 
on L, the eigenvalues of the two L-dependent electronic states (solid circles) change as L 
changes. 
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